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SUMMARY 

An approximate thermal equation was derived for quenching distance 
hased on a previously proposed diffusional treatment. The quenching 
distance is expressed in terms of the thermal conductivity^ the fuel 
mole fraction^ the heat capacity^ the rate of the rate -controlling chem- 
ical reactions^ a constant that depends on the geometry of the quenching 
surface^ and one en^irical constant. 

In order to use the equation with experimental data, the rate- 
controlling reaction must he specified. Two choices were made: (l) the 

reaction between active particles and fuel (corresponding to the diffu- 
sional equation), and (2) the reaction between oxygen and fuel. 

The effect of pressure on quenching distance was shown to be in- 
versely proportional to the pressure dependence of the primary flame 
reaction, with a small correction necessitated by the effect of pressure 
on flame temperature . 

The general equation was used with the Semenov equation for burning 
velocity to show that the quenching distance was inversely proportional 
to burning velocity and pressure at any given initial temper at Tn:e and 
equivalence ratio. 

The two equations, one based on the active-particle - fuel reaction, 
the other based on the oxygen-fuel reaction, were tested, using pub- 
lished data for the quenching of propane -oxygen-nitrogen flames. The 
data selected included the effect on quenching distance of oxygen- 
nitrogen ratio, propane concentration, and initial mixture tenrperature 
and press-ure. A correlation of these data was obtained, using each of 
the two equations, although both possessed shortcomings. The equation 
using the active -particle - fuel reaction did not correlate data for all 
rich mixtures; the equation using the oxygen-fuel reaction, while cor- 
relating data for both rich and lean mixtures, showed a larger deviation 
from the predicted linear relation than the other equation. 
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INTRODUCTION 

The process of flame quenching is of interest^ since it may he re- 
lated to other combustion phenomena of engineering iirportance, such as 
flame stabilization, flammability limits, and the general behavior of 
flames near cold walls. Flame -quenching processes become especially im- 
portant in turbojet combustion systems when operation at low pressure 
(i.e. , high altitudes) is considered. Flame quenching is usually studied 
experimentally in terms of the quenching distance, which is defined as 
the minimum channel size that will allow a given flame to propagate. 

Most of the available quenching data have been obtained either by (l) 
observation of the minimum tube diameter or rectangular slot width that 
will allow a flame to flash back, or by (2) determination of the minimum 
distance between plane -parallel plates that will allow a flame to propa- 
gate from a spark of minimum ignition energy. 

A rigorous theoretical treatment of the quenching process appears 
to be a most difficult task (refs. 1 and 2). Consequently, present 
theoretical treatments of quenching are necessarily approximate and 
seek primarily to correlate the quenching process with the variables 
that affect it, such as press-ure, temperature, fuel type and concentra- 
tion, inert-gas concentration, and quenching-STirface geometry. 

Approximate treatments of quenching have been based either on a 
thermal or a diffusional mechanism for the process (refs. 3 and 4), 

Such apparently different outlooks yield results consistent with experi- 
mental data because the equations for heat and mass transfer are for- 
mally identical, and because the thermodiffusivity and the molecular 
diffusivity are numerically nearly equal (ref. 5). Thus, each purely 
thermal approach to a combustion process has a diffusional analog. 

A quenching-distance equation that successfully correlates lean 
quenching data (refs. 4, 6, and 7) is the diffusional equation given by 
Simon, Belles, and Spakowski (ref. 4). Unfortunately, the use of this 
equation in its present form is limited to stoichiometric or lean mix- 
tures. The source of this difficulty appears to be the choice cf the 
reaction kinetics. 

Essentially, the objective of this paper is to extend, if possible, 
the useful range of the quenching concepts of reference 4 to include 
hydrocarbon-rich mixtures. Because of the similarities between heat 
and mass flow, it is to be expected that a thermal analog equation can 
be derived that will correlate data equally as well as the diffusional 
equation of reference 4. It was believed that such an alternate thermal 
equation might be more susceptible to changes in the reaction kinetics 
than the original diffusional model. It was hoped that the result of 
such changes would be an extension of the useful range of the equation 
to include rich mixtures. 
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This report contains the derivation of a thermal quenching equation 
in which the rate -controlling reaction is not specified. Tsvo possible 
rate -controlling reactions are post\fLated, and the resulting two quench- 
ing equations are tested using published quenching data. It is shown 
that, if the oxygen-fuel reaction is assumed to be rate -controlling, 
quenching data are satisfactorily correlated for both rich and lean 
propane -oxygen-nitrogen flames. The pressure dependence of the quench- 
ing distance and the relation of quenching distance to burning velocity 
are also discussed. 


THECRY 

As an initial assun 5 >tion, the authors of reference 4 assumed that 
in order for a flame to propagate the nirnber of reaction events per 
cubic centimeter that occur in the gas ahead of the bTorning zone (here- 
inafter called the reaction zone) must be above some minimum value. 

The number of reaction events per cubic centimeter in the reaction zone 
is assumed to be the nimiber of active particles per cubic centimeter 
Cj^ m\altiplied by the average number of effective collisions (the 

average active-particle chain length) made by each particle before its 
destruction at the wall. The minimum number of reaction events neces- 
sary for flame propagation is assumed to be a constant fraction A of 
the total number of molecules per cubic centimeter Nip. In terms of 

these definitions, the criterion for quenching is given as 

VpCp < ANp (l) 

for a single species of active particle and a surface chain-treaking 
efficiency of 1 (symbols are defined in the appendix). 

The definition of the number of reaction events necessary for flame 
propagation used in reference 4 is not a unique one. Another possibility 
is to assume that the necessary number of reaction events in the reaction 
zone is a constant fraction F of the total number of events per cubic 
centimeter x that normally occur during the passage of a flame through 
a unit volume of gas. The quenching criterion then becomes 

< Fx (2) 

This quenching criterion leads to a more suitable thermal equation than 
would be obtained through the use of equation (l). 

Equation (2) may be converted to its thermal counterpart by multi- 
plying through by hRT/P to give as a final result (expressed in 
cal/mole) : 
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Hr < FHt (3) 

Here^ the number of calories per mole produced in the reaction zone 

must be equal to or less than some critical value, which is assumed to 
be a constant fraction F of the total heat produced per mole of un- 
burned gas Bp. 

For the combustion of hydrocarbons, complete combustion of the fuel 
is assumed to occur for all mixture compositions. Consequently, the 
total heat produced by the process may be written as the product of the 
mole fraction of fuel in the unburned gas Xp and the heat produced 

upon the disappearance of 1 mole of fuel AH (as defined herein, AH is 
constant and equal to the heat of combustion to COg and H2O only for 
stoichiometric or lean mixtures j for rich mixtures, AH decreases as the 
oxygen concentration decreases). Thus, for hydrocarbon flames, equation 
(3) may be written as 


Ep<FXpAH (4) 

The expression given in reference 4 for the chain length is 


V , = 


1 GiDiT. 


(5) 


The thermal egui valent of this equation is found by multiplying through 
by hCj^RT/P to obtain as a final result 

d^hc.RT 


= 


G^D.t.P 

ill 


According to the definitions given in reference 4, 


( 6 ) 



(7) 


The heat produced per reaction event in the reaction zone is proportional 
to the heat produced upon the disappearance of one fuel molecule. For 
simplicity, the two are assumed to be equal so that the heat produced 
per reaction event may be written as 


h = 


N 


( 8 ) 
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Since the molecular diffusivity and thermodiffusivity are asstimed 
to be equivalent. 


D-, = 


x^RT 


PC. 




(9) 


Any small error in 5 >licit in the assumption of equation (9) should not 
influence the over-all behavior of the resulting equation for quenching 
distance . 


Combination of equations (6) to (9) yields the thermal analog of 
equation (S) : 


Hr = 


d^ AH C 




V 


G^Nx^. 


( 10 ) 


Simtiltaneous solution of equations (4) and (8) gives for the thermal 
quenching equation 


d 


2 


0p,r 'W 


( 11 ) 


(The inequality sign has been removed since the quenching distance is 
related to the tube configuration that just quenches a given flame.) The 
most ing)ortant fact concerning this equation is that the form for w is 
no longer specified; any reaction may now be chosen as rate -controlling. 


If equation (l) rather than equation (2) had been used as the 
quenching criterion, the final thermal quenching equation would have 
been formally identical with equation (ll) , except that the fiiel mole 
fraction Xj. would not appear. This difference is unin^jortant as re- 
gards the correlation of data if only a narrow range of fuel concentra- 
tion is considered. For data that cover a wide range of fuel concen- 
tration, the thermal quenching equation as written in equation (ll) is 
found to be the more satisfactory of the two possibilities. 


In terms of the active-particle - fuel reaction (the reaction spec- 
ified as rate -controlling in the diffusional equation of ref. 4), the 
reaction rate w is given by 


- kt^-r ^i,r 
i 


w 


( 12 ) 
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If a colllsion-theory-type temperatOTe dependence is assigned to 
equation (l2) becomes 




(13) 


Combination of equations (lO) and (l3) gives 

^2 ^ ^ F 


C B a V c exp(.5t) 

Lj i.r r / 




(14) 


p^r t r i^r r 


The choice of rate -controlling reaction is not limited to the 
active -particle - fuel reaction. For example^ it is possible to follow 
Semenov (ref. 8) and assume the reaction in the combustion zone to be 
bimolecular and first-order with respect to fuel and oxygen. For this 
assumption^ the real reacting species are obviously not fuel and oxygen 
molecules^ but such an assumption is satisfactory if the concentrations 
of the reacting species are proportional to the fuel and oxygen con- 
centrations. In this case^ then, the rate of the reaction is given by 

w == (l5) 


If a collision-theory-type temperature dependence is assigned to 
equation (l5) becomes 


1/2 

w = Bga^b^,!^/ expU 


£ 
.RTr 


(16) 


The combination of equations (lO) and (l6) gives 


d^ = 


FG.Nx X„ 
X r f 






/-E, 


r-^r 


exp 


li 




F 1 

Bo ^s 


(17) 


In order to test equations (l4) and (l7), data for the quenching 
of propane -oxygen-nitrogen flames were used. These data include the 
effect on quenching distance of oxygen-nitrogen ratio, fuel concentra- 
tion and pressure (ref. 6), and initial temperati^re (ref. 9). 

In order to use equation (l4), the assuirptions made in reference 4 
concerning the active -particle - fuel reaction were followed: T^. was 

taken to be 0.7 T^, the partial pressure of fuel in the reaction zone 
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was assumed to be one -half the partial pressure of fuel in the unburned 
§as^ the partial pressures of the active particle-s (H atoms, OH radi- 
cals, and 0 atoms) in the reaction zone were assumed to be 0.7 times 
their partial pressure in the adiabatic equilibrium flame and was 

assumed to be 7 kilocalories per mole (ref. lO) , 

In order to use equation (l7), the partial pressures of fuel and 
oxygen in the reaction zone were taken as equal to their values in the 
unburned gas and Eg was assumed to be 38 kilocalories per mole (ref. 

ll) . In keeping with previous usage (refs. 4 and 6), was chosen 
to be eqml to 0.7 Tp». 

The thermal conductivities were calculated for the unburned gas 
mixture at 0.7 Tp. The conductivities for oxygen, nitrogen, and propane 

were calculated at 1553^ K, using the tables and data given in chapter 
8 of reference 12. Conductivities at other temperatures were computed 
on the assumption of a 3/4-power temperature dependence. The thermal 
conductivities of the gas mixtures were calculated as the sum of the 
mole fraction times the thermal conductivity of each component, or, 

X = Xj^Xi (18) 

i 

It was felt that this simple linear mixing rule was adequate for systems 
largely composed of nitrogen and oxygen because of their similar 
conductivities. 

The heat capacities were calciolated for the unburned gas mixture 
at 0.7 Tp by the same linear mixing rule used for the thermal con- 
ductivities. Individual heat capacities for oxygen, nitrogen, and pro- 
pane were calculated from data given in reference 13. 

The equilibrium adiabatic flame temperat-qres and product composi- 
tions were calculated by the matrix method of reference 14 using the 
thermodynamic constants of reference 14 and the heat of formation of 
propane given in reference 13. 


RESULTS Am DISCUSSION 

Correlation of Quenching -Distance Data 

Inasmuch as the square of the quenching distance is predicted to 
be directly proportional to (eq. (l4)) or to (eq. (l7)), values 

of 'ilf^ and were calculated for the pressure, equilibrium adiabatic 

flame composition and temperature, unb-urned gas composition and temper- 
ature, and tube geometry (assumed to be plane -parallel plates), which 
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correspond to the quenching distances reported in references 6 and 9 
for propane -oxygen-nitrogen flames. Plots’ of d^ against and 

are shown in figures l(a) and (b), respectively. The ranges of the var- 
iables included were: pressure, 0.1 to 1.0 atmosphere; equivalence ratio, 

0.33 to 1.90; mole fraction oxygen in the oxidant mixture, 0.21 to 0.70; 
and unburned gas temperature, 300° to 558° K. In figure l(a), it is 
seen that (which is entirely analogous to the diffusional equation 
of ref. 4 except for the choice of quenching criterion) satisfactorily 
correlates the effect of these variables except for propane-rich mix- 
tures where (p > 1.2. 

Inspection of figure l(b) reveals that ijfg (which includes the 

oxygen-fuel reaction as the rate -determining step, rather than the 
active -particle - f\ael reaction) correlates not only the propane-lean, 
but also all the propane-rich quenching data. Thus, it appears that by 
proper choice of the rate-controlling reaction, the quenching concepts 
of reference 4 may be extended in such a way as to include rich mixtures 
for the propane -oxygen -nitrogen system. 

A mean line drawn thro\igh the data presented in figure 1 reveals 
that the predicted linear relation between d^ and and ilfg does 

not hold exactly. Rather, the data are best correlated if 

d^ oc (■'If. )*^'^^ oc (\jf )0*®'^. This difficulty is related, at least in part, 

"G S 

to the choice of reaction mechanism^ since the "best correlation” ex- 
ponent is changed by a change in reaction mechanism. 


Pressure Dependence of Quenching Distance 

Examination of the equation for quenching distance (eq. (ll)) re- 
v*eals that all terms on the right side are independent of pressure ex- 
cept the reaction rate w. However^ hydrocarbon flame tenperatures 
change slightly vith pressure because of changes in the extent of dis- 
sociation. Consequently^ a pressure change indirectly affects 
temperature- and coirposition-dependent terms in the equation. This 
effect is quite small^ so that the pressure dependence of the quenching 
distance may be regarded as primarily dependent on the pressure de- 
pendence of the initial chemical processes in the combustion vave. 

The calculated pressure dependence of the quenching distance (in- 
cluding the indirect effect of pressure on flame temperature and compo- 
sition) is conpared in table I with the experimental values of reference 
6. The average deviation of the calculated values from the observed 
values was found to be 10 percent for and 17 percent for \jrg. 
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Reference 6 gives a value of 7 percent for the average deviation of 
values computed for the diffusional equation. 

It is seen that ijfg does not predict the press\ire dependence as 
well as does except for the high-velocity flames at an oxygen 

fraction of 0.70. This is probably because the assiaraption of a simple 
second-order reaction between fviel and oxygen is a poor one. Concern- 
ing prediction of pressixre dependence, a better choice of over-all order 
wovild be 1.7 rather than 2.0; this is in agreement with theoretical 
studies of the pressure dependence of burning velocity (ref. 12, p. 765), 
which indicate that the global reaction in most hydrocarbon flames ranges 
between first and second order. Such a choice would also tend to improve 

the linearity of the relation between d^ (observed) and \|r (calculated) 

o 


Quencliing Distance and Bnrning Velocity 


The relation of quencliing distance to burning velocity has been 
discussed in references 3, 4, and 15. In reference 4^ an equation re- 
lating the two variables is derived on a purely diffusional basis. The 
thermal analog to this equation may be derived as follows : 


According to reference 16^ the Semenov equation for the burning 
velocity may be written as 






(19) 


Multiplying equation (ll) by equation (l9) and writing for a^ and 
Po^ 


WX^P 


(20a) 


Po 



(20b) 


yields 


RT 


d = 


UP 




/^l\^/xRT - 

2x x_ FG.I — I W 

r F i\n 2 j 




C C w 
p,F 


( 21 ) 
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If it is assumed that W and w have the same pressure dependence, 
it follows from equation (21) that, at a given equivalence ratio, 

U oc ioc (22) 

From reference 6, it is seen that n decreases from about -0.9 to 
-1.05 as the percent oxygen in the oxidant mixture is increased from 21 
to 70. Simultaneoiisly, the burning velocity increases from about 40 
centimeters per second to the order of 300 centimeters per second. 
Consequently, equation (22) predicts that the escponent describing the 
pressure dependence of U shoixLd Increase from about -0.1 to 0.05 as 
U increases from 40 centimeters per second to about 300 centimeters 
per second. This result may be conpared with the experimental findings 
of Lewis (ref. 17) (which indicate an increase of about -0.05 to 0.07 
as burning velocity changes from 40 to 300 cm/sec) to confirm equation 
(22) qualitatively. 


CONCLUDIITG REMARKS 

The fact that reasonable agreement exists between experiment and 
quenching concepts based on either heat or mass transfer Indicates only 
that these processes have the same formal laws and that the thermal and 
molecular diffusivities are either equal or remain directly proportional 
for the data considered in this study. No evidence for the preponder- 
ance of either heat or mass transfer in quenching can be deduced from 
these results. The principal advantage gained by formulation of the 
quenching concepts of reference 4 in a "thermal" form is that a change 
in reaction kinetics is much simpler from a conceptual point of view 
for the "thermal" form than for the "diffusional" form. This con- 
ceptual advantage might be useful in a discussion of bizarre fuel- 
oxidant systems for example, CS2-F2 flames. In addition, the actual 

calculations may be simplified, since by the proper choice of reaction 
mechanism, it may become unnecessary to calculate the equilibrium 
flame con^josition. 


SUMMARY CF RESULTS 

The results of an investigation into the possibility of extending 
quenching concepts previously proposed to include quenching of 
hydrocarbon-rich flames may be summarized as follows: 

1. An equation for the critical quenching configuration of a 
channel was derived for hydrocarbon flames: 
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where 

Cp^j, heat capacity in reaction zone, cal/ (°K)( mole) 

d characteristic dimension of tube geometry; quenching distance, 

cm 

F constant that relates total number of reaction events that occxir 

during the passage of a flame to number that must occ\ir in the 
reaction zone for flame propagation 

dimensionless geometrical factor, dependent only on channel 
geometry 

N Avogadr o ' s nimiber 

w rate of reaction in reaction zone, molecules/ (cc) (sec) 

Xp mole fraction of fuel in unbiirned gas 

mean thermal conductivity in reaction zone, cal/(cm) (sec) (°K) 

2. Two possible flame -initiating reactions were considered in the 

detailed formulation of w: (l), the reaction between active particles 

and fuel molecules, and (2), the reaction between oxygen and fuel 
molecules . 

3. The two equations res\ilting from the two reaction mechanism 
choices were tested using published data that included the effect of 
oxygen -nitrogen ratio, fuel concentration, pressure, and unburned gas 
tenqjerat-ure on the quenching distance for propane -oxygen-nitrogen flames 
These variables were correlated reasonably well for both rich and lean 
flames by the equation involving the oxygen-fuel reaction. The equation 
involving the active-particle - fuel reaction was satisfactory primarily 
for lean flames. 

4. It was concluded that, by the proper choice of reaction kinetics 
the quenching concepts previously proposed may be cast in a form useful 
for the prediction of both propane-rich and -lean quenching data. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, December 6, 1954 
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APPENDIX - SYMBOLS 

The following symbols are used in this report: 

fraction of molecules present that must react for flame to con- 
tinue to propagate 

fuel concentration^ molecules/cc 

Arrhenius constant 

oxygen concentration^ molecules/cc 

concentration of i^^ active particle^ molecules/cc 

heat capacity, cal/ (^K) (mole) 

average heat capacity, to cal/ (^K) (mole) 

diffusion coefficient, cm^/sec 

characteristic dimension of tube geometry; quenching distance, 
cm 

activation energy, cal/mole 

constant that relates the total number of reaction events that 
occur during the passage of a flame to the number which must 
occur in the reaction zone for flame propagation 

dimensionless factor, value of which is dependent only on tube 
geometry 

heat released upon consumption of 1 mole of fuel by combustion 
process, cal/mole 

heat produced by chemical reaction in primary reaction zone , 
cal/mole 

total heat produced by combustion of 1 mole of unburned gas , 
cal/ mole 

heat produced per reaction event, cal 
rate constant, (cc) (molecules)/sec 
average molecular weight of unburned gas, g/mole 
molecularity of flame reaction 


•zcn Q 
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N 


% 

n 


P 

R 

T 

U 

W 


Xf 


X 


a 

X 

V 

p 

T 

<P 




ATOgadro’s number 

total number of molecules present per cc 

exponent describing pressure dependence of quenching distance 

moles of reactant per moles of product from stoichiometric 
eqmtion 

pressure^ atm 

gas constant 

temperature y 

flame speedy cm/sec 

average reaction rate in flame front as defined by Semenov, 
molecules/ (cc) (sec) 

reaction rate in primary reaction zone, molecules/(cc) (sec) 
mole fraction of fuel 

number of reaction events per cc that occur during passage of 
flame through a unit volume 

mole fraction of oxidant in oxidant-inert mixture 

thermal conductivity, cal/(cm) (sec) (^K) 

average active-particle chain length 

density, g/cc 

time between effective collisions, sec 
equivalence ratio 
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Subscripts : 

F flame 

i active -particle species 

o unburned gas 

r primary reaction zone 

s involves reaction of oxygen and fuel molecules 

t involves reaction of active particles and fuel molecxiles 
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TABLE I. - COMPARISON OF OBSERVED PRESSURE DEPENDENCE 
OF QUENCHING DISTANCE WITH PREDICTED VALUES 


Oxygen 

fraction^ 

a 

Equivalence 

ratio^ 

<P 

Exponent n describing 
pressure dependence, 

d CX p-n 

Observed^ 

Predicted 

Diffu- 

sional^ 

Thermal 



0.17 

0.943 

0.90 

0.88 

0.83 

1.03 


1.000 

.89 

.85 

.76 

1.01 


1.340 

.84 

.75 

.75 

1.00 


1.530 

.71 

.75 

.75 

1.00 

0.21 

0.738 

0.85 

0.86 

0.80 

1.00 


.864 

.84 

.89 

,84 

1,03 


1.000 

.89 

.93 

.89 

1.03 


1.240 

.95 

.78 

.79 

1.02 


1.490 

.98 

.76 

.77 

1.00 

0.30 

0.566 

0.76 

0.87 

0.87 

1.01 


.662 

.93 

.89 

.88 

1.04 


.773 

1.06 

.92 

.90 

1.07 


1.000 

.98 

.94 

.90 

1.09 


1.380 

.93 

.86 

.84 

1.06 


1.903 

.74 

.76 

.76 

1.01 

0.50 

0.476 

1.01 

0.91 

0.89 

1.06 


.544 

1.01 
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Figure 1. - Summary correlations of quenching data by thermal quenching equations (data 
from refs. 6 and 9). 
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Figure 1. - Concluded. Summary correlations of quenching data by thermal quenching equations 
(data from refs. 6 and 9). 
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